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Noisy metamolecule: strong narrowing of fluorescence line 
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We consider metamolecule consisting of bosonic mode correlated with the two-level system: it 
can be, for example, plasmonic mode interacting with the quantum dot. We focus on the parameter 
range where all the correlations are strong and of the same order: interaction between bosonic mode 
correlated with the two-level system, external coherent drive and dissipation. Quantum Monte-Carlo 
simulations show a fluorescence of this system at dissipation larger than the driving amplitude and 
strong (by the order of magnitude) narrowing of its spectral line. This effect may be related to kind 
of a quantum stochastic resonance. We show that the fluorescence corresponds to finite domain over 
the coherent drive with sharp lower threshold and there is splitting of the Wigner function. 


Correlated bosonic modes and discrete spectrum sys¬ 
tems attract much attention last time: they are key el¬ 
ements of applications in quantum optics chem¬ 

istry [ 6 ], and biology [7] that includes nanolasers [8HI3], 
qubits, quantum shuttles [14], cold atoms [T5|, plasmonic 
nanoprobes and biosensors. m Today technology en¬ 
ables to reduce the size of the bosonic mode host up to 
nanometer scale. In particular, such plasmonic nanosys¬ 
tems are often referred to as “metamolecules” HI]- 
[ 22 |. Here we consider as metamolecule also plasmonic 
nanosystems interacting with quantum dot [23] . 

In most applications, correlated bosonic modes and 
discrete spectrum systems in the first approximation 
may be reduced to the Jaynes-Cummings class of mod¬ 
els. They are quite well investigated except the param¬ 
eter range where there is no well defined small or large 
interaction parameter or there is high level of dissipa¬ 
tion and there are large fluctuations. After Ref. [24| 
we refer to these range as to the “no-man’s-land” [25] . 
In these cases advanced computational efforts are re¬ 
quired na w\- For applications mentioned above this 
parameter range is quite typical. In the real systems the 
noise is always present. Usually this is the problem espe¬ 
cially for applications related to quantum effects. Recent 
experiments on single DNA hairpins show amplification 
of the response by stochastic noise [28l |29]. Here we 
show partially motivated by experiments that noise may 
effectively increase coherence in metamolecule while the 
system occupies the “no-man’s-land”. 

In the simplest case the Jaynes-Cummings model 
(JCM) corresponds to the noisy bosonic mode correlated 
with the two-level system [30l|3T]. Electromagnetic mode 
in the resonator or the plasmon in the nanostructure 
may serve as the bosonic mode while ground and ex¬ 
cited states of a molecule or a quantum dot may serve 
as the two-level system m m 0 Ea. In the “no-man’s- 
land” the Rabi frequency is of the order of dissipation in 
bosonic mode. From the fluctuation-dissipation theorem 
dissipation means noise. We do nonequilibrium quan¬ 
tum Monte-Carlo simulations [261 EH E] of the system 
and show that under coherent drive there is a kind of a 


dynamical “phase transition” driven by noise accompa¬ 
nied by by the narrowing of the fluorescence spectral line 
by the order of magnitude, enhancement of the signal 
to noise ratio and splitting of the Wigner function — it 
plays the role of the effective “order parameter”. From 
experimental point of view this fluorescence regime may 
find wide applications where the strong noise is the un¬ 
avoidable problem. 

We consider within JCM metamolecule, see Fig.[2 The 
dynamics of the system in Markov approximation is gov¬ 
erned by the Lindblad equation [34] : 


dp/dt 


—i 



1- ^ (2dpd^ — pa^d — d^dp) + 

y {2apa+ - pa+a - a+ap) , (1) 


where 7 ^ and 75 are the dissipation rates of the two- 
level system and bosonic mode correspondingly and H = 
H}) T Hq^ T Va—w T ^b—w T ^a—b' Here Hi) — uJiycHd is 
the Hamiltonian of the bosonic mode (d is annihilation 
operator), Ha = is a two-level atom Hamiltonian 

(d = 1 ^) (e| is operator of transition from the excited 
state |e) to the ground one |^)), I^-a = U^(d+(j + d+d) 
is the interaction between the bosonic mode and the two- 
level system (TLS) where CLr is the Rabi frequency. The 
coherent drive of TLS and bosonic mode is described by 
Va-w = + d+) and V^-w = ^b{d~^ + d). We choose 

the system of units here and below where h = 1. 



FIG. 1. (Color online) Possible realisation of noisy meta- 
molecule consisting of plasmon in metallic nanoparticle laiisi 
na (or nanorod 13 E]) correlated with the two-level system 
(atom, molecule or a quantum dot). 
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Eq. 0 has been intensively investigated last decades 
both analytically and numerically. However this is not so 
for the “no-mans-land”, where any analytical progress is 
very hard. (In [26] “no-man’s-land” has been investigated 
but with slightly different model and parameter range.) 
In addition, in this parameter region the system is in be¬ 
tween quantum and classical regimes since the average 
boson number, N = (d'^'d), is larger than unity but not 
big enough for applicability of the semiclassical methods 
or 1/N expansion [35]. Note that experimentally meta¬ 
molecule is realized in nanolaser type systems. 

In Ref. [ 8 ] the plasmonic nanolaser has been experi¬ 
mentally realised. This experiment is realisation of JCM 
where the role of the bosonic mode plays the plasmon 
in the metallic nanoparticle, see Fig. while TLS cor¬ 
respond to quantum dots. For this system 75 ^ 10^^ — 
Due to strong field localization in the plasmonic 
mode the Rabi interaction constant may become com¬ 
parable with the damping in the plasmonic mode [36] . 
Qr We should note that in this ex¬ 

periment many TLS have been correlated with plasmonic 
bosonic mode and in our model — only one. The other 
difference is the external drive: in the experiment the 
system was driven by 5ns pulses while we consider co¬ 
herent drive. We do not see however any problem with 
the realisation of coherent drive by an external laser and 
the use of a single quantum dot. The realistic ampli¬ 
tude of the laser field is ^ 10^ — lO^R/m. For realistic 
semiconductor quantum dots typical relaxation rate is 
7a ^ 10^ — 10^^ 5 “^ |8]. Another possibility is quantum 
dot or molecule in the cavity [DiaElISll like we sketch in 
Fig. [l}i. More recent realization corresponds to biology 
where coherent effects in the presence of large noise have 
been observed m 

So, we imply the following values of constants: 7 ^ = 
and Qr = The dumping rate of 

bosonic mode 75 G (0.5-which corre¬ 
sponds to different nanocavity materials. Also we choose 
O 5 = lOQa = ^ which corresponds to more intensive in¬ 
teraction between external field with nanocavity. It takes 
place when plasmonic mode is excited, e. g., see Ref. jS] 



FIG. 2. (Color online) “Lasing” threshold, a) boson number 
A, b) Mandel parameter (red solid) and coherence function 
(black dash). Maximum corresponds to the threshold where 
narrow line fluorescence appears. Insert shows in arbitrary 
units the signal to noise ratio (SNR). It is inversely propor¬ 
tional to the boson mode linewidth. 
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FIG. 3. (Color online) “Lasing” threshold and above, (a) -(c) 
The steady-state reduced bosonic density matrix in the basic 
of the occupation numbers. Its multimodality corresponds to 
the narrow line fluorescence regime (“lasing”), (d) -( f) show 
the Wigner function where x and y are real and imaginary 
part of the coherent state. Break up of the Wigner function 
into two peaks corresponds to narrow line fluorescence regime. 
Graphs (a),(d); (b),(e); and (c),(f) were prepared at Q = 
0.04, 0.06, 0.08. Density plots (bootom raw) show the Wigner 
function at (g) Q = 0.03, (h) Q = 0.05 (In the region at the 
center of the “attol” the Wigner function is negative, i.e. this 
is non-classical state), (j) D = 0.07 and (k) Q = 0.09. 


The absolute value of interaction constant Q. with exter¬ 
nal field we will vary from 5 • to 1 • which 

corresponds to experimentally achievable laser field am¬ 
plitude, 10 — lOOR/m. In simulation we normalize 
all values on the So we consider the “no-mans- 

land” case when 7 ^ VLr ^ D ^ 7 ^. Also we suppose 
the exact resonance between bosonic mode and two-level 
system. 

What happens beyond these parameters? In the insert 
of Fig. Dd we show the signal to noise ratio; in fact, any 
point within its broad extremum ridge corresponds to 
narrow line fluorescence. 

To solve Eq. 0 we use the quantum Monte Carlo 
algorithm where the Lindblad equation reduces to the 
piecewise deterministic process on quantum trajecto¬ 
ries [26l |33l [34] , with the Gilbert space dimension N ^ 
100; this is more than enough, see Fig. I^a)-(c), where 
the peaks show the most probable N. 

First we investigate the evolution of the average boson 
number N with the drive amplitude D to detect the pa¬ 
rameter range where N grows that might correspond to 
narrow line fluorescence. We plot N in Fig. [^ not only 
in the quantum case, corresponding to the exact solution 
of the master equation 0 . but for comparison we show 
in Fig. 0 results of well-known approximate solutions 
for A: 1) the semi-classical approximation and 2) “non- 
interracting” case corresponding to the dynamics of the 
bosonic mode in the external field but without interac¬ 
tion with the TLS. 

The semiclassical solution of Eq. 0 implies total ne¬ 
glecting of quantum fluctuations. As follows from Eig. 0 
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FIG. 4. (Color online) Spectra evolution of the bosonic mode: Q = 0.05,0.06,0.07,0.08,0.102,0.104,0.106,0.108,0.12,0.13. 


there is hysteresis (bistability) in this case. Exact so¬ 
lution (red curve) of 0 shows no hysteresis: it disap¬ 
pears due to strong quantum fluctuations [24]. Out¬ 
side the bistability the value of population inversion, 
([d+,(j]_) > 0. So TLS can be pumped using just 
a monochromatic external wave and dissipative bosonic 
mode in agreement with the “bad-cavity limit” [38| . 

In our problem quantum fluctuations play important 
role as shows the Mandel parameter Q = — 

{N)‘^)/{N) — 1 = {N){g^‘^\0) — 1), where N = d+d and 
g^‘^\r) = (d+(t + r)d+(t + r)d(t)d(t)). Q has maximum 
at the driving field amplitude corresponding to the be¬ 
ginning of narrow line fluorescence regime, see Fig. It- 

Fig-i shows linear growth of N above certain thresh¬ 
old of drive. Below we investigate the parameter range of 
the transition. In Figs. (a)-(c) we present the reduced 
density matrix in the occupation number representation 
for different values of the coherent field amplitude Q. At 
small f], much below the transition to narrow line fluores¬ 
cence, there is one pronounced peak in the density matrix 
diagonal corresponding to = 0. At the transition and 
above the density matrix becomes essentially nondiago¬ 
nal. We see the appearance of the phase multimodality 
of the matrix elements in the narrow line fluorescence 
regime: i.e., splitting of the density matrix distribution. 
In Figs. (a)-(c) the density matrix splits at least on 
three ridges. This behaviour of the density matrix im¬ 
plies essentially nonthermal state of the system. Below 
we shall probe how coherent are bosonic modes and ir¬ 
radiation in this state. At larger values of ft the density 
matrix again becomes diagonal. This is the end of the 
narrow line fluorescence regime. 

More informative is the Wigner function evolution pre¬ 
sented in Figs.[^(d)-(f) for the same driving amplitudes. 
At very small driving fields the Wigner function has 
only one peak located near the state with zero quanta 
of bosonic mode. Above the threshhold, Fig. i (h), it 
has the “atoll” form and the blue hole corresponds to 
negative values of the Wigner function. This is known 
for dissipative JCM under coherent excitation of bosonic 
mode. m But here we consider the generalised Hamil¬ 
tonian where also coherent driving of TLS is taken into 
account. Then we see new behavior: the Wigner function 
first splits up into two narrow peaks (narrow line fluores¬ 
cence regime), Fig. [^j), but going to larger field ampli¬ 


tudes we observe that the left peak of the Wigner func¬ 
tion disappears and only the right one survives, Fig. 
(k). This behaviour should correspond to the end of nar¬ 
row line fluorescence regime. 

Finally we consider the spectra of the bosonic correla¬ 
tors, S{t) = (d+(t + r)d(t)). Graphs in Fig. [^explicitly 
illustrate the presence of the lower and upper threshold 
of the domain where the linewidth do reduces. 

Generally speaking there are a number of known cases 
when noise may effectively increase coherence: this is the 
stochastic resonance PlliO]. Metamolecule is strongly 
correlated system and it is not clear in advance what 
noise may do. Possible interpretation of the narrowing 
of fluorescence line is the quantum stochastic resonance. 
However it can be a kind of lasing as well: We have the 
threshold and the population inversion in our system. We 
see that spectral line, which length is determined by the 
level of spontaneous emission, is narrowing. This points 
out that coherence the system output signal increases. 
Another point is splitting and narrowing of the Wigner 
function. The characteristic width of the Wigner func¬ 
tion determines quantum fluctuations. In the regime of 
the developed generation the same situation takes place. 
In Ref. Ell three lasing criteria have been formulated: 1) 
transition to lasing should be accompanied by nonlinear 
dependence of the boson number on external drive, 2) 
emission linewidth should strongly decrease in the lasing 
regime and 3) the second order coherence function g^‘^^ (0) 
should tend to one which indicates the coherent output. 
In our model fluorescence in the regime of “narrow spec¬ 
tral line” shows all these properties. 

Also we note that the metamolecule we focus on has 
many features typical for recent nanolaser type experi¬ 
ments where the bosonic mode is the plasmonic mode of 
metallic nanoparticle while semiconductor quantum dot 
or dye molecule play the role of the active medium and 
similar narrowing of spectral line and threshold behav¬ 
ior have been seen laiiiii]. In these experiments pulse 
pumping has been used. However the usual duration of 
the pump pulse was about 10 ns. While all characteristic 
times of nanolaser/spaser like problems are much shorter 
(the nanosecond scale). So this pulse pumping can be 
well approximated by coherent pumping like we do. Re¬ 
cently another experiment showed the fluorescence under 
coherent pumping identified with lasing 02 ]. 
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To conclude, we describe the behavior of metamolcule 
driven by the external field. It is shown that fluorescence 
at dissipation larger than driving is accompanied by nar¬ 
rowing of the spectral line by the order of magnitude and 
splitting of the Wigner function of the emitted light. This 
effect may find applications in quantum systems where 
strong noise due to high losses is unavoidable problem 
like in plasmonics. 
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